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ABSTRACT 

We investigate the evolution of supermassive black hole mass (Mbh) and the host 
spheroid mass (M sp h) in order to track the history of the MBH-M sp h relationship. 
The typical mass increase of Mbh is calculated by a continuity equation and accretion 
history, which is estimated from the active galactic nucleus (AGN) luminosity function. 
The increase in M 8p % is also calculated by using a continuity equation and a star 
formation model, which uses observational data for the formation rate and stellar 
mass function. We find that the black hole to spheroid mass ratio is expected to 
be substantially unchanged since z ~ 1.2 for high mass objects (Mbh > 10 8 - 5 M Q 
and M sp h > 1O 113 M0). In the same redshift range, the spheroid mass is found to 
increase more rapidly than the black hole mass if M sp h > 10 11 M Q . The proposed 
mass-dependent model is consistent with the current available observational data in 
the MBH-M sp h diagram. 
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1 INTRODUCTION 

Recent observations show that most, if not all, nearby massive spheroids (elliptical galaxies, lenticular and spiral bulges) have 
supermassive black holes(SMBHs) with masses in the range of Mbh = 1O 6 -1O 1O M0 at their center. I n the local universe, the 
SMBH mass is correlated wi t h characteristi c parameters of the host spheroid: the luminosity L.,„h jKormendv fe Richstonel 



19951 ; iMarconi fe Hunt] 120031 ; iGrahaml l2007h, the mass M SDh of the spheroid jMagorrian et al ll99Sl ; Irlaring fe RixIboO^ ') 



the stellar veloc ity dispersion ff dGebhardt et aljf 



20od: iFerrarese fe MerrlfruboOC 



Trcruamc et al. 2002) and the various other 



properties (e.g.. iGraham et al.l (|2007l ) : I Aller fc Richstonel (j2007l h iKisaka. Koiima fc Otanil (|2008l )h Understanding the origin 
of these relations is thought to provide insights into the "co-evolution" of SMBHs and the host spheroids. The role of central 
black holes in galaxy formation and evolution is not yet established since the SMBH mass is tiny compared with the galaxy 
as a whole. Areas tha t are understood at present, as well as those suggested for future observation, were recently reviewed in 
Cattaneo et al. ( 20091 ) . Some tentative theories have been proposed to clarify the origin of these correlations (e.g., Silk fc Reesl 



(1998 



and they are fairly successful at reproducing the correlations at the present time. However, no theoretical model is 
widely accepted since there are many unknown parameters in the evolutionary model. For exa mple, in semianalytic models, 
depen ding on the model assumptions, the evolution of the scaling relation differs markedly (e.g., Croton ( 20061 ) ; Hopkins et al 
(2003)). Therefore, it would be of value to study the evolution using an alt ernative appro ach 



Cosmic evolution of the ratio Mbh / M sp h is observationally studied. McLure et al. (2006) used the 3C RR sample of 
radio-loud active galactic nuclei (AGNs) for < z < 2. The number of sources for each redshift bin is not sufficiently large to 
examine the 2-dependence. By fitting data for the relation Mbh I M sv h oc (1 + z) 1 , the index is estimated as 7 = 2, although 
data for z < 1 are also consistent with no evolution (7 = 0). iTreu et al.l (2007) also derived a similar result for SMBH masses 
in a sample of 20 Seyfert galaxies around z = 0.36; however, the evolution was slightly weak, Mbh I M gv h oc (1 + z ) 1 ' 5 . It is 



not easy at present to judge whether or not cosmological evolution exists, due to several selection biases (jLauer et al.l l2007) 



* E-mail: kisaka@theo.phys.sci.hiroshima-u.ac.jp 
f E-mail: kojima@theo.phys.sci.hiroshima-u.ac.jp 



2 S. Kisaka and Y. Kojima 



The evolution of the correlations can also be inferred from the demographics of galaxies and AGNs 1 Merloni. Rudnick fe Di Matted 



2004). In some reasonable theories, the nuclear activity of AGNs is proportional to the mass accretion r ate into the S MBHs. 



The luminosity function of AGNs and the redshift evolution measure the buildup of SMBH m ass (e .g., ISoltanl (jl982j)). The 
evolution of the galaxies is related to the star formation history. Merloni. Rudnick fe Di Matted (2004) combined the evolution 
of the SMBH mass density derived from the mass accretion history with that of the stellar mass density derived from the 
star formation history of the Universe, in order to investigate the evolution of the MBH-M ap h relation. They found that the 
growth of SMBHs appears to predate that of spheroids as Mbh / M sp h oc (1 + z) 0,4-0,8 , for which the power index of the 
redshift factor is smaller than that obtained by observations. 

Here, we consider the evolutionary model by taking into account the mass-dependent growth in the phenomenological 
analysis. That is, a simple power-law relation M bh / M sv h oc (1 + z) 7 , is not a s sume d. Only integrated quantities have 
been used so far in phenomenological studies (e.g., Merloni. Rudnick fc Di Matted ( 2004 )). Our model is an improvement of 
previous models, in which the mass-dependent effect is now included as an important factor. The mass-dependent property 
of cosmological evolution is known as " d ownsi zing", in which star formation becomes active in less massive galaxies as the 
cosmic time increases fe.g JCowie et all (|l996h ). Similarly, an increase in the number of AGNs with redshift compared with 
the local universe is followed by a decline be yond a peak whose r edshift depends o n luminosity and he nce on the SMBH mass, 
assuming a constant Eddington ratio (e.g.. lYencho et al. Jiool)). More recently. Habita et all (j2009h have studied the mass 
dependence of redshift for active black holes by using the Malmquist-bias-unaffected method, and argued that AGN samples 
at lower redshifts are increasingly dominated by less massive black holes. These observations suggest that the evolution of 
both SMBH mass and galaxy mass depends on their masses. In this paper, we extend the evolutionary models based on the 
AGN luminosity function and star formation history by including the mass-dependent effect and derive the evolution of the 
MBH-M sv h relation at certain cosmic ages. 

We organ i ze the paper as follows. In Section 2, we br i efly r eview the evolutionary model of SMBHs considered by 
Marconi et all |2004h : Ishankar. Weinberg fc Miralda-Escudel (|2009h . in order to clarify the model parameters. We assume a 



mean Eddington ratio for the active SMBHs. By adopting a continuity equation and an assumed initial condition, the AGN 
luminosity function correlates directly with the SMBH mass function at all times. In Section 3, we describe the evolution 
model of spheroids. Using a method similar to that used for SMBHs, that is, a continuity equation and an assumed initial 
condition, the stellar mass function and specific star formation rate are used to calculate the mass increase of the spheroids. In 
Section 4, we combine these two evolutionary models and discuss the evolution of the MBH-M ap h relation. Finally in Section 
5, we provide a summary. Throughout the paper, we adopt a cosmology with f2 m = 0.3, Qa = 0.7 and Ho = 70 km s _1 
Mpc" 1 . 



EVOLUTIONARY MODEL OF BLACK HOLES 



In the following section, we review the evolu tionary equation of the SMBH mass function according to lMarconi et al.l (|2004l ); 
Shankar. Weinberg fc Miralda-Escudel ( 20091 ). The mass of black holes generally increases through mass accretion, which also 
causes AGN activity. The SMBH m a ss function is theref o re det e rmined from the obse r ved luminosity function of AGNs (e.g., 
Cavaliere. Morrison fc Wood! l|l97ll ): ISmall fc Blandfordl |l992h : lMarconi et all l|2004h : IShankar. Weinberg fc Miralda-Escude 
|2009h l. We assume that AGNs are powered by accretion into SMBHs, and that SMBH growth takes place during a phase 
in which the AGN is shining at a fraction A of the Eddington luminosity, converting the mass accretion rate with a radiative 
efficiency e. The number density of SMBHs with mass Mbh in co- moving space is denoted by u(Mbh ,t), which satisfies the 
continuity equation 



dn(M BH ,t) 
dt 



+ 



t) 



8Mb 



[ti(Mbh , t) (Mbh {Mbh , t))] = 0, 



(1) 



where {Mbh{Mbh ,t)) represents the mean mass-growth rate of Mbh at time t. The mass growth i s related to AGN activity 
and hence is expressed by the bolometric luminosity function of the AGNs. Marconi et al. 1 20041 ) derived an evolutionary 
equation of the SMBH number as 



dn(M BH ,t) 
dt 



(1 - e)\ 2 c 2 



d$(L bo i,t) 



dLboi 



(2) 



where tsdd is the Eddington time. 

The mass growth of SMBHs between redshifts z\ and 22 is calculated by the integration 



M(z 2 ) - M( Zl ) = 



{M BH (MBH,t))^dz 



(1 - e)A 
rffidd In 10 



-rrr-t [$(L 6o ; , t)]^dz. 

n(MBH,t) dz 



(3) 
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The evolutions of SMBH mass density and typical mass at z are calculated by using eqs. ([2]) and ([3]), which are 
expressed by the AGN luminosity function with two parameters e and A. We adopt the AGN l u minos ity function of 
Shankar, Weinberg fe Miralda-Escudi |2009l ) with a smoothing modification by iRaimundo fc Fabian! (2009). The observed 



AGN luminosity function is converted to the bolometric luminosity by using the luminosity-dependent correction o flMarconi et al 
' 20041). The initial condition for the integr ation of eq. (3) is chosen as the condition at z = 6, which is the same as 



Shankar. Weinberg fc Miralda-Escudi (|2009l ). They noted that the integration does not significantly depend on the initial 



value after z ~ 3.5 unless the SMBH duty cycle is extremely small. 

In order to constrain the free parameters e and A, we compare the local SM BH mass function derived from the evolution 
model with an observational value. Shankar. Weinberg fc Miralda-Escudi ( 20091 ) carried out a detailed comparison by using the 
SMBH mass function derived by various methods. We also obtained the same results as lShankar. Weinberg fc Miralda-Escude 
(2009), even when including the technical smoothing modification. The best-fitting values are e = 0.065 and A = 0.42. 
Uncertainty in the radiative efficiency e is not a highly important parameter affecting the mass growth of SMBHs since it 
affects only the overall normalization. The parameter is fixed as e = 0.065 from now on. The Eddington ratio A is a critical 
parameter describing SMBH evolution. Shankar. Weinberg fc Miralda-Escudi 1 20091 ) estimated a reasonable range of their 
model parameter, which is converted as A = 0.28-0.56. The range implies not uncertainty in any statistical sense, but rather 
a preferable range. 

In Figure Q3 we plot the averaged growth history of SMBHs with different masses for < z < 1.2. The mass range of 
black holes for which current observational data are available is limited to 10 7 < Mbh /Mq < 10 9 ' 5 . It is clear that the growth 
depends on black hole mass. A black hole with Mbh ~ 1O 7 M0 rapidly increases its mass from z — 1.2 to z = 0, while one 
with Mbh ~ 10 Mq stops its growth at z = 1.2. We also demonstrate the upper and lower limits of the Eddington ratio 
for two representative evolutionary tracks as shaded areas. One is the track of Mbh = 1O 8,5 M0 at z = 0, and the other is 



uncertainty does not seriously affect the evolution for Mbh > 10 s Mq. The uncertainty in the track of less massive black holes 
can be clearly seen, for example, ~ 0.4dex in log Mbh for Mbh = W 7 ' 5 Mq at z = 0.9. 



EVOLUTIONARY MODEL OF SPHEROIDS 



In this section, we introduce a simple model of spheroid mass evolution. It is well known that galaxies form two sequences in 
color-magnitude space: star-forming, late- type galaxi e s occu py the blue cloud, whereas quiescent, bulge-dominated early-type 
galaxies reside on the red sequence (e.g.. iBell et al.l (|2004l )). We focus on early type galaxies and derive the spheroid mass 
evolution assuming a fixed bulge-to-t otal ratio (B/T). 

We use the method reported by IBell et all (|2007l ) to determine the evolution of the stellar mass function. The number 
densities of both the early- and late-type galaxies, respectively, satisfy the continuity equation with a source function S for 
each type as 



dn earl y(M,t) + [n ear i v (M, t)M ea Hy {M, t)] = S early (M,t), 
dn latE {M,t) + _^_ [nZote(Mit)Mtate(Mjt)] = Slate (M,t). 
For all galaxies, we have 

^"jf^ + -^[n all {M,t)M aU {M,t)] = S all {M,t), 
where 

n a u(M,t) = (M,t) + n late (M, t), 

n a ii(M, t)M a ii(M, t) — it-early 

(M, t)M eaH y(M, t) + m a te(M, t)M late {M, t), 
Sall(M,t) = Searly(M, t) + Sl a te(M,t). 



(4) 
(5) 

(6) 

(7) 
(8) 
(9) 



Bell et al.l (|2007l ) calculated the evolution of the mass function for all galaxies with an assumption S a ii = 0, and found that 



the results were consistent with the observational COMBO-17 data at 0.2 < z < 1.0 (jBorch et al]|2006l) . They also found 
dni a te/dt — by comparing the observed mass function of all galaxies with that of late-type galaxies between z — and 
z ~ 0.9. We adopt both conditions in our analysis. The continuity can be written as 

dn ear i y (M ,t) d 
dt ~ ~dM l 

d_ 

dM l 

where 



[n a „(M, t)M a it(M, t)] 

[nearly {M ear ly{M,t))], 



(10) 
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(M earlv (M,t)) = J^iMJL M all{M ,t). 

n ea rly(M,t) 

The typical mass of the early-type galaxies is obtained by integrating the equation 



M ear ly(Z2) — M ear ly(Zl) 



{M ear i y {M,t))^-dz. 

dz 



(11) 



(12) 



The spheroid mass evolution is determined from th e stellar mass function of ear ly-type galaxies adopting B/T=0.7, which 
does not depend on the redshift or galaxy mass fsee lTamura. Ohta fc Uedal |2006h ). 



The mass increase M a ii is determined by the star formation rate. We adopt the staged r model introduced bv lNoeske et al 



I 2007) and use the same values and uncertainty ranges. The main mode l parameters are c a and cp , which control the gas 
exhaustion timescale and the formation redshift, respectively. The model bv lNoeske et al.l (|2007l ) is consistent with observations 
over a wide range of masses and may be used as a well-fitting formulation of specific star formation rate (SSFR). Figure [2] 
sh ows the SSFR ^ /M at z =0.3, 0.5, 0.7 and 0.9. Here we include much more observational data, and obtain the same results 



Noeske et al. 1 2007 ). For the stellar mass function data, we use the data of Ilbert et al. ( 2010h from the COSMOS survey 



covering 2-deg , in which morphological and spectral classifications were carried out. The combined classification s allow us to 



isolate the "blue elliptical" galaxies that are not included in the red sequence. We use the early-type (elliptical in lllbert et al 
1 20101 )1 galaxy mass function at z < 1.2 as the spheroid mass function. 



The mass evolution is determined by easily integrating eq. (|12[) since eq. (|lip is expressed as a smooth curve. The results 
are shown in Figure[3] It is found that there is no significant growth in the mass range of M > 1O 11 M0. As discussed regarding 
the evolution of SMBHs, uncertainty in the model parameters is important in the phase of rapid growth. The parameters of 
spheroid evolution ar e the B/T ratio , and the SSFR parameters (c a and cp). We use B/T=0.7, but examine the variation for 
0.4< B/T <1, which llm et al.l (|2002f ) considered as a reasonable range. The range of variation is shown for two evolutionary 
tracks in Figure[3] The uncertainty in log M sp h for M sp h = 1O 1O ' 3 M0 at z = 0.9 is ~ 0.6dex. This uncertainty can be neglected 
for M ap h > 1O 11 M0. We also examine the variation due to the SSFR parameters in Figure [4] It is found that the variation in 
the SSFR parameters affects the evolution more strongly than that in the B/T ratio. However, neither affects the spheroid 
evolution for M aph > 10 11,5 M Q . 



RESULTS 



4.1 Evolution of the MBH-M ap h relation 

The combination of models for SMBH mass and spheroid mass allows us to investigate the mass -dependent evolution of the 
MBH-M ap h relation. In the local universe, the Mbh-M^ relation found bv lHaxing fc Rbd ([2004]) satisfies 



log(M M /M ) = 8.20 + 1.121og(M sph /10 11 M Q ). 



(13) 



Using this, we consider the relation at < z < 1.2. Figure [5] shows our results for a reference model with A — .42, B/T=0.7, 



10- 



In order to account for the intrinsic scatter at z = 0, we add an offset of ±0.3dex (H aring fc Rix 



c a = W M -' and cp 

2004) in eq. (|13p . so that two curves corresponding to upper and lower limits are also plotted at each value of z. We found 



that the MsH-M ap h relation does not change since z ~ 1.1 in the region for Mbh ^ 10 Mq and M ap h ^ 10 M©. The 
reason is clear in Figures [T] and O where no evolution was seen in either mass on the massive side. On the other hand, there is 
a significant deviation from the linear relation when Mbh < 1O 8 M0. This denotes that a massive BH with Mbh = 1O 8 M0 is 
already located at a relatively small spheroid mass of M aph = 10 10 - 5 M Q at z ~ 1. The spheroid mass subsequently increases 
to M ap h = 1O 1O ' 8 M0, while Mbh is fixed at 1O 8 M0. In order to evaluate the increase of the SMBH mass for a fixed spheroid 
mass M ap h, we consider the difference Alog(Af_B£r) = \og(MBH(z) j Mbh(0))- When the evolution is plotted as log(l + z) 7 , 
we obtain 7 = 2.11 at M aph = 10 10 M Q , 7 = 1.05 at M aph = 10 10 ' 5 Af Q and 7 = 0.43 at M aph = 10 n M Q . In the higher-mass 
regime of M aph > 10 n M Q , no significant evolution is found. Thus, the evolution of the MsH-M ap h relation at z ^ 1.1 clearly 
depends on the mass range. The growth is more rapid for a smaller spheroid mass. If we integrate the relation with the mass 
function over the whole mass range, then the averaged power index becomes 7 ~ 0.8. We show the relation with 7 = 0.8 
for comparison in subsection 4 . 2 (Fig ure El). This parameterization for evaluating cosmological evolution has been used by 
Merloni. Rudnick fc Pi Matted (2004), and our result is consistent with their v alue within an error of la. The i ndex of the 
average is almost the same as that of the previous phenomenological model (e.g.. iMerloni. Rudnick fc Di Matted (|2004r 0. but 
mass-dependent evolution, especially at lower mass, is a remarkable new property in our model. 

In order to compare the growth of SMBHs and spheroids, we present two extreme cases in Figure [6] One is calculated 
without evolution of SMBHs, and the other is calculated without evolution of spheroids. Both two curves deviate from our 
reference model as well as from the MBH-M sp h relation at z — 0, as z increases. The deviation is remarkable at lower mass. 
However, our reference model is close to the curve for no evolution of SMBHs. This means that the mass increase of the 



spheroids around M. 



ph 



10 Mq is more rapid. 
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We examine the model parameters discussed in Section 2 and 3 to demonstrate deviation from the reference model. 
The three parameters are independent and the results change monotonically as each parameter is changed. Thus, we have 
calculated two extreme cases: upper and lower limits of the MBB-M sp h relation. These results are shown in Figure [Jj For 
the model with c a — 10 20 ' 7 and cp — 10 _3 '°, which corresponds to the curve for the upper limit in Figure [Jj the distribution 
of Mbh in the range of M sp h 10 11 Mq is a nearly horizontal line at z ~ 1.1. On the other hand, for c a — 10 20 ' 4 and 
cp — 10 -1 ' 7 , corresponding to the curve for the lower limit in Figure [jj the distribution in the same range is even lower than 
that of the local relation. For lower masses, the deviation from the linear line is general, and the typical SMBH mass depends 
on the model. For example, our model shows that Mbh ranges from 10 7 M Q to 1O S - 5 M for M sph = 10 10 ' 5 M Q at z = 1. The 
possible range of masses is wide, but our reference model, which can be regarded as an average of two extreme cases, suggests 
slightly larger SMBH mass than that of the local universe relation. It is also found that the relation is not affected in the 
mass range of Al sp h ~ 10 Mq. 



4.2 Comparison with observations 

In this subsection, we compare our model with observational results available in the literature. McLure et al. ( 20061 ) investi- 
gated the evolution of the SMBH-to-host galaxy mass ratio in the redshift range of < z < 2. They estimated the masses 
Mbh for the 3C RR quasars and spheroid masses M av h for the 3C RR radio galaxies. Two samples are not the same, but the 
averaged mass was determined. Recently, Decarli et al. ( 2009bh studied the redshift dependence of the relation up to z — 3 
using a sample of 96 quasars for which the host galaxy luminosity is known. They listed Mbh and M sp h for each source. Both 
results are used for comparison after we derive Mbh from the observati o nal data by the s ame m ethod. Black hole masses 
are usually estimated by the virial method (e.g., McLure fc Jarvis ( 20021 ) ; Kollmeier et al. ( 2006l )V The masses are derived 
from the broad emission-line width v as a velocity indicator, and the monochromatic continuum luminosity L is used as an 
indicator of the region size: 



log Mbh = A + 2 log v + B log L, 



(14) 



where A and B are constants. The line width v is typically taken from the broad line of Mgu or H/3, and the luminosity L 
is taken at 3000Aor 5100A. The constants A and B have not been determined with high accuracy, since the best fit values 
slightly depe nd on observa t ional choic es of v and L . The bl ack hole masses are derived by a different choice of A and B in 
the works of iMcLure et all (j2006h and becarli et all (|2009bl ). It is prefera ble to use a single formula for the estimation. One 
possible approach is to use the values, A and B, derived by McGill et al. ( 20081 ') . who cross-calibrated a number of different 
formulae in a redshift range where more than one broad emission line could be observed simultaneously in optical spectra. 
We adopt their calibrated val u es for A a nd B, in order to re-d erive black hole masses. We calculate Mbh of the objects in 
the samples by McLure et al. (200(3) an d Decarli et al. ( 2009bh . for which the line width of Mgu or H/3 is measured. The 
corresponding line width of nine objects is not listed; consequently, these are excluded from our analysis. In Figure we 
plot (Mbh, M S ph) not for each sour ce bu t for the average, sinc e we focus only on statistical properties. Our values Mbh and 
those listed in M cLure et al. |2006h and iDecarli et al l (|2009bl are not notably different by averaging. iBennert et all (|201oh 
investigated the e volution of the MBH -Lsvh relation. □ As discussed previously, Mbh is derived from spectral data using the 



mass estimator bv lMcGil 



et al 



( 2008h. The s phero id mass M sp h is derived from L ap h and the mass-to-light ratio by a certain 
stellar population model! Bruzual fc Chariot (|2003h ). Thus, converted data are also plotted in Figure [8] We group data into 
two ranges, M sp u < lO n M0 and M sp u > 1O 11 M0, and respectively average Mbh and M sp h, since the sample number is 
small. Although detailed comparison is difficult due to the small sample and the ambiguity in the conversion relation, the 
results are consistent with mass-dependent evolution. The samples are deawn from the high mass end of the spheroid mass 
function (M sp h > 1O 11,5 M0), where we expect little or no evolution in the MBH-M sp h relation since z ~ 1. Comparing our 
evolutionary model with observations in each redshift bin we find consistency with observations, although we cannot exclude 
a mass-independent scenario because of the lack of observational data in the low spheroid ma ss range. 



We next consider additional observation data, although the comparison becomes indirect. ISalviander et al.l (|2007r ) inves- 
tigated the evolution of the Mbh-o relation. An evolutionary model of a is necessary, but here we assume that the conversion 



the relation between a and M sp h at z — is derived 
and the relation is assumed to 



relation from a to M 3p h is the same as that of z = 0. In other words, 

by eliminating Mbh in eq. (|13(l and the Mbh-o relation at z = bv lTremaine et al. 

hold for all z. Averaged data points from lSalviander et al. (2007) are also plotted in Figure(8] These data correspond to small 
masses of M sp h < lO n M0 and show relatively strong evolution, supporting the hypothesis of a mass-dependent evolution 
scenario. 



Similar samples were previously given bv lTreu et al. | J2007h , but most of these were covered by the sample of Bennert et all j2010l ). 



6 S. Kisaka and Y. Kojima 



5 SUMMARY AND CONCLUSIONS 

We have modeled the mass-dependent evolution of black holes and spheroids since z ~ 1.2. The black hole mass evolution 
was investigated by using a continuity equation and the observed AGN luminosity function. The spheroid mass evolution was 
also derived from a continuity equation and the observed SSFR. Both evolutions, which were consistent with "downsizing", 
significantly depend on the masses of the SMBHs and spheroids, and the Mbh-M 3P h relation is mass-dependent. We assumed 
that all galaxies at z = satisfy the linear relation of the log(MBff)-log(M sp h) diagram with some uncertainty, and examined 
the history to z ~ 1.2. It was found that the relation is unchanged for < z < 1.2 in the range of Mbh > 1O 8 M0 and 
M sp h > 10 11 Mq. In the low-mass regime, however, there is clear deviation from the linear relation. The mass increase of the 
spheroids is likely larger than that of black holes, although there is uncertainty in the model parameters. Our model suggests 
that SMBHs were already located in spheroids with relatively small mass, compared with the local relation. 

We also compared our results with currently available observations concerning the relations between Mbh and host 
spheroid parameters at certain cosmic ages. Our mass-dependent evolutionary model is consistent with these observations, 
although the observational sample is not sufficiently complete to make a statistically significant comparison. Larger obser- 
vational samples are needed to study the history of SMBHs and their hosts statistically, expecially in the low-mass regime 
where the differences between the mass-dependent and the mass-independent evolution models are higher. If a break from 
the linear log(Maiy)-log(M sp ft) is found, our mass-dependent scenario would be supported and it would be possible to put 
stronger constraints on model parameters. It would also be of great interest to extend the present results to higher redshift, 
where the model uncertainties are considerably smaller. This would enable us to reveal possible inconsistencies between model 
and observations. However, the extension of the present study to hi gher redshift is no t straightforward, as the adopted stellar 
mass function for early-type galaxies is only significant at z < 1.2 l llbert et al. 2O10l ). 



REFERENCES 

Aller, M. C, & Richstone, D. O. 2007, ApJ, 665, 120 

Bell, E. F., Zheng, X. Z., Papovich, O, Borch, A., Wolf, O, & Meisenheimer, K. 2007, ApJ, 663, 834 
Bell, E. F., et al. 2004, ApJ, 608, 752 
Bennert, V. N., et al. 2010, ApJ, 708, 1507 

Borch, A., Meisenheimer, K., Bell, E. F., Rix, H.-W., Wolf, O, Dye, S., Kleinheinrich, M., & Kovacs, Z. 2006, A&A, 453, 
869 

Bruzual, G., & Chariot, S. 2003, MNRAS, 344, 1000 

Buat, V., et al. 2008, A&A, 483, 107 

Cattaneo, A., et al. 2009, Nature, 460, 213 

Cavaliere A., Morrison, P., & Wood, K. 1971, ApJ, 170, 223 

Chen, Y.-M., Wild, V., Kauffmann, G., Blaizot, J., Davis, M., Noeske, K., Wang, J.-M., & Willmer, C. 2009, MNRAS, 393, 
406 

Cowie, L. L., Songaila, A., Hu, E. M., & Cohen, J. G. 1996, AJ, 112, 839 
Croton, D. J. 2006, MNRAS, 369, 1808 

Damen, M., Labbe, L, Franx, M., van Dokkun, P. G., Taylor, E. N., & Gawiser, E. J. 2009, ApJ, 690, 937 

Decarli, R., Falomo, R., Treves, A., Labita, M., Kotilainen, J.K., & Scarpa, R. 2009, MNRAS, accepted (arXiv:0911.2988) 

Ferrarese, L., & Merritt, D. 2000, ApJ, 539, L9 

Gebhardt, K., et al. 2000, ApJ, 539, L13 

Graham, A. W. 2007, MNRAS, 379, 711 

Graham, A. W., Driver, S. P., Allen, P. D., & Liske, J. 2007, MNRAS, 378, 198 
Haring ,N.,& Rix, H.-W. 2004, ApJ, 604, L89 

Hopkins, P. F., Hernquist, L., Cox, T. J., Keres, D. & Wuyts, S. 2009, ApJ, 691, 1424 
Iglesias-Paramo, J., et al. 2007, ApJ, 670, 279 
Ilbert, O., et al. 2010, ApJ, 709, 644 
Im, M., et al. 2002, ApJ, 571, 136 

Kisaka, S., Kojima, Y., & Otani, Y. 2008, MNRAS, 390, 814 
Kollmeier, J. A., et al. 2006, ApJ, 648, 128 
Kormendy, J., & Richstone, D. 1995, ARA&A, 33, 581 

Labita, M., Decarli, R., Treves, A., & Falomo, R. 2009, MNRAS, 396, 1537 
Lauer, T. R., Tremaine, S., Richstone, D., & Faber, S. M. 2007, 670, 249 
Magorrian, J., et al. 1998, AJ, 115, 2285 
Marconi, A., & Hunt, L. K. 2003, ApJ, 589, L21 

Marconi, A., Risaliti, O, Gilli, R., Hunt, L. K., Maiolino, R., & Salvati, M. 2004, MNRAS, 351, 169 



Evolution of the M BH -M sph relation 7 



Martin, D. C, et al. 2007, ApJS, 173, 415 

McGill, K. L., Woo, J.-H., Treu, T., & Malkan, M. A. 2008, ApJ, 673, 703 
McLure, R. J., & Jarvis, M. J. 2002, MNRAS, 337, 109 

McLure, R. J., Jarvis, M. J., Targett, T. A., Dunlop, J. S., & Best, P. N., 2006, MNRAS, 368, 1395 
Merloni, A., Rudnick, G., & Di Matteo, T. 2004, MNRAS, 354, L37 
Noeske, K. G., et al. 2007, ApJ, 660, L47 

Pannella, M., Hopp, U., Saglia, R.P., Bender, R., Drory, N., Salvato, M., Gabasch, A., & Feulner, G. 2006, ApJ, 639, LI 
Perez-Gonzalez, P. G., et al. 2008, ApJ, 675, 234 
Raimundo, S. I., & Fabian, A. C. 2009, MNRAS, 396, 1217 

Salviander, S., Shields, G. A., Gebhardt, K., & Bonning, E. W. 2007, ApJ, 662, 131 

Shankar, F., Weinberg, D. H., & Miralda-Escude, J. 2009, ApJ, 690, 20 

Silk, J., & Rees, M., J. 1998, A&A, 331, LI 

Small, T. A., & Blandford, R. D. 1992, MNRAS, 259, 725 

Soltan, A. 1982, MNRAS, 200, 115 

Tamura, N., Ohta, K., & Ueda, Y. 2006, MNRAS, 365, 134 
Tremaine, S., et al. 2002, ApJ, 574, 740 

Treu, T., Woo, J.-H., Malkan, M. A., & Blandford, R. D. 2007, ApJ, 667, 117 
Ueda, Y., Akiyama, M, Ohta, K., & Miyaji, T. 2003, ApJ, 598, 886 
Yencho, B., Barger, A. J., Trouille, L. & Winter, L. M. 2009, ApJ, 698, 380 

Zheng, X. Z., Bell, E. F., Papovich, C, Wolf, C, Meisenheimer, K., Rix, H.-W., Rieke, G. H., & Somerville, R. 2007, ApJ, 
661, L41 



8 S. Kisaka and Y. Kojima 




Figure 1. Averaged growth history of SMBHs computed using AGN luminosity function and e = 0.065. Eddington ratio is A =0.42. Red 
shaded areas show the difference in logMgu between A = 0.56 and A = 0.28 for representative tracks at z =0.3, 0.5, 0.7, 0.9 and 1.1. 
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Figure 2. Specific star formation rate as a function of stellar mass from z = 0.3 (upper left), z = 0.5 (upper right), z = 0.7 (lower left) 
and z = 0.9 (lower ri ght). Thick lines sh ow staged r models and thin- lines show effect o f varying c a and Ca (or, zj and t ). Fitting pa- 
ramet ers are given in lNoesk e et al.l j2007l). Data poi nts correspond t o measurements fro m lzheng et al.l J2007l)(open c ircles ), iMartin et al.l 
l|2007h (clo sed circles)T]l glesias- P^amo" et al.l l|2007l ) (open squares). iBuat et al.l l|200Sf l (closed squ ares) , |Bell e t al.l ll2007l )(open upward 
tria ngles') , iPannella et al.l l|2006j) (closed upward triangles) , iDamen et al. I j2008l ) (open diamonds), IChen et al.l ^2009^) (closed diamonds) 
and|P erez-Gonz^Iez^TTu |2008h (open downward triangles) . 
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Figure 3. Averaged growth history of spheroids computed using B/T=0.7 and SSFR with c a = 10 20 ' 7 and cp = lCT 2 ' 7 . Red shaded 
areas show the difference in logM sp ^ between B/T=0.4 and B/T=L0 for representative tracks at z =0.3, 0.5, 0.7, 0.9 and 1.1. 
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Figure 5. Evolution of Msn-^sph relation from z ~ 1.2 to z = 0. Model parameters are A = 0.42, B/T=0.7, c a = 10 20 ' 7 and 
cp = 10 -2 ' 7 . An offset ±0.3dex is added due to scatter in the local Mgfj-Mgph relation. Relation, and upper and lower limits at z = 
are plotted by thick dotted lines and thin dotted lines, respectively. 
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Figure 6. Comparison of evolutionary models of MsH-^sph relation. Black solid line is our reference model; black dotted line is 
relation at z = 0. Model with fixed SMBH mass is denoted by top (red) line; model with fixed spheroid mass is denoted by bottom (blue) 
line. 
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Figure 7. Evolution of MsH-Mgph relation from z ~ 1.2 to z = varying model parameters. Black solid line is our reference model; 
black dotted line is M BH -M sph relation at z = 0. Top (red line) is upper limit, for which A = 0.28, B/T=0.4, c a = 10 20 ' 4 and cp = 10 -1 ' 7 ; 
Bottom (blue line) is lower limit, for which A = 0.56, B/T=1.0, c Q = 10 20 7 and cp = 10 -3 . 
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Figure 8. Comparison of MBH-^sph relation from z ~ 1.2 to z = with observations. Solid lines arc referen ce model s and dotted lines 
are local relations. The re l ation shifted by (1 + z)° 8 is shown by dash-d otted (green) lin es. Red circles: I Bennert et all ((2010*); light blue 
squares: ISalviander et al. | J2007D : purple triangles: [becarli et all i2009bh : blue diamonds: iMcLure et al. | j2006h . Error bars are standard 
deviation of the logarithmic mean. 



